Global entanglement in XX Z chains 
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We examine the thermal entanglement of XXZ-type Heisenberg chains in the presence of a 
uniform magnetic field along the z axes, through the evaluation of the negativity associated with 
bipartitions of the whole system and subsystems. Limit temperatures for non-zero global negativities 
are shown to depend on the asymmetry A but not on the uniform field, and can be much higher 
than those limiting pairwise entanglement. It is also shown that global bipartite entanglement may 
exist for T > even for A > 1, i.e., when the system is fully aligned (and hence separable) at T = 0, 
and that the bipartition leading to the highest limit temperature depends on A. 
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I. INTRODUCTION 

Quantum entanglement is a fundamental trait of quan- 
tum mechanics, representing the non-local correlations 
with no classical analogue that can be exhibited by com- 
posite quantum systems. Interest on entanglement has 
been triggered in the last few years by the discovery of its 
potential for developing new forms of information trans- 
mission and processing [lHj|, being a resource in Quan- 
tum Information science [.5| . Entanglement has also pro- 
vided a new perspective for the analysis of correlations 
and transitions in many-body quantum systems (6l-[l0|. 
In systems at finite temperature T, two fundamental 
questions which immediately arise [ill, 113 ^^^i a) the 
determination of the limit temperatures for the existence 
of different kinds of entanglement and b) the possible 
emergence of entanglement for T > when the ground 
state is separable due to entangled excited states. 

Let us first recall that a mixed state p of a bipartite 
quantum system A + B is said to be separable or classi- 
cally correlated if it can be expressed as a convex combi- 
nation of product densities [l^, i.e., p — QoiPa ® Pb^ 
with qa > 0, J2a 9ct = 1 and p^, p% density matrices for 
each component. Otherwise, p is entangled or insepara- 
ble. Separable states can be created by local operations 
and classical communication (LOCC) and satisfy Bell in- 
equalities as well as other classical properties such as be- 
ing more disordered globally than locally [13, [3 ■ Pure 
states (p^ = p) are separable if and only if p = pA® Pb^ 
but this is not necessarily the case for mixed states, where 
the determination of separability is in general an NP hard 
problem [l^, [l3| • Accordingly, while for pure states bi- 
partite entanglement can be measured in terms of the 
entropy of the reduced density of a subsystem |18|, a 
rigorous computable measure of entanglement for mixed 
states has so far been obtained just for the case of two 
qubits, through the concurrence [19]. 

Nonetheless, bounds for entanglement in general mixed 
states can be obtained by means of the negativity [20I - I23 , 
which is a measure of the degree of violation of the crite- 
rion of positive partial transpose (PPT) 0,[23| in entan- 
gled states and is easily computable. Although the PPT 



criterion is a necessary separability condition, sufficient 
just for two-qubit or qubit-|-qutrit systems, the negativity 
fulfills some fundamental properties of an entanglement 
measure [20*1 , being an entanglement monotone and pro- 
viding bounds to the teleportation capacity and distil- 
lation rate. In n-qubit systems at finite temperature, it 
can then be employed [23 to detect the entanglement be- 
tween the components of any bipartition {m} — {n — to} 
of the whole system, as well as of any subsystem, be- 
yond the level of pairwise two-qubit entanglement mea- 
sured by the concurrence. Another fundamental result 
is that any state p is completely separable (convex com- 
bination of n-product densities) if it is sufficiently close 
to the fully mixed state I/d [2l|, [H, [l^l {d is the sys- 
tem dimension) . This ensures that any canonical thermal 
state p{T) cx exp[— i7/T] of a finite system described by 
a Hamiltonian H becomes separable above a finite limit 
temperature, since it will be as closed as desired to I/d 
for sufficiently high T . In particular, no bipartite entan- 
glement can arise if Tr(p - I/df < [d{d - 1)]"! 

The aim of this work is then to examine, by means 
of the negativity, the global bipartite entanglement of 
thermal states p{T) of 71-spin chains interacting through 
an XXZ-iype coupling placed in a transverse magnetic 
field. Heisenberg chains can be employed for solid state 
quantum computers [28^1 and XXZ chains have been 
used to describe quantum computers based on NMR [2^ 
and on electrons on Helium [sO, [Ml . Relevant studies of 
the pairwise thermal entanglement between two qubits in 
Heisenberg chains have been made [i[l3,,29, 32 36], re- 
vealing rich phenomena. In small XXZ chains [29|,l3J] it 
has been shown that there is no pairwise entanglement for 
T > for anisotropics A > 1, i.e., when the ground state 
is fully aligned for any field, and that the corresponding 
limit temperature remains bounded for large negative A 
if n is odd (vanishing V A in the n = 3 antiferromagnetic 
case at zero field 34]). In contrast, we will show here 
that global entanglement between m and n — to qubits 
may also exist for A > 1 if n > 4 and T > 0, implying a 
thermal reentry. Moreover, limit temperatures for non- 
zero global negativities do not saturate for A — >■ —00 in 
odd systems, but diverge in both the ferro- and antiferro- 
magnetic cases, even for n = 3. These temperatures are 
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independent of the magnetic field but dependent on A 
and the bipartition, displaying crossings as A is varied. 
Numerical results up to n = 10 qubits will be provided, 
together with full analytical results for n — 3. 



II. FORMALISM 

We will consider an XX Z Hamiltonian for n qubits or 
spins of the form 

H ^ bS,- jy^is^si + + v^slsil (1) 

where denotes the spin operator at site i, Sz — ^1 
is the total spin z-component, b accounts for the Zeeman 
coupling to a uniform magnetic field and s^s^ ~^ ^y'^y ~ 
(s!|_siL + s!_si!|_)/2 is the hopping or entangling term. Our 
attention will be centered on a cyclic chain with nearest 
neighbor coupling (w^ = VaSj^i+i for a — x, z, with n + 
1 = 1), although some results for the fully connected case 
(f Jj* = Vo,y i < j) will also be commented for comparison. 
In the first case the spectrum of H and the entanglement 
of its eigenstates are independent of the sign of Vx for even 
n, as it can be changed by a rotation around the z axes 
at odd sites {s^^y — >■ (— l)*s^^y). In any case, [H, Sz] — 0, 
so that the eigenstates of H have good quantum number 
M (eigenvalues of Sz)- 

We will consider the n-spin thermal state 

p{T) = Z-'cM-H/T], (2) 

where Z — Tr exp[—H/T] is the partition function and 
T > the temperature (we set Boltzmann constant k — 
1). Global entanglement between m and n — m selected 
qubits will be analyzed by means of the negativity [13] 

iV,[p] = iTr(|p*^|-l), (3) 

where p = {m} — {n~ m} denotes the bipartition and p*'' 
the ensuing partial transpose of p [1^. Eq. (jS]) is just the 
absolute value of the sum of the negative eigenvalues of 
/9*p (as Trp'j" = 1), so that according to the PPT criterion 
[H,!!!], Np[p] > indicates entanglement between the 
m selected qubits and the rest. Eq. satisfies as well 
properties of an entanglement measure [2Q|] , being a con- 
vex function of p and an entanglement monotone (it does 
not increase under LOGO). Moreover, distillable entan- 
gled states satisfy Np[p] > [l^. Although the present 
analysis leaves away bound PPT entangled states and 
does not capture all aspects of n-partite entanglement 
(separability of all bipartitions does not necessarily im- 
ply full separability), it goes beyond the standard analy- 
sis based on pairwise two-qubit entanglement (the latter 
may vanish even if Np[p]> for all global bipartitions, 
as occurs in GHZ states [231). For pure states p = |^')(\E'|, 
it can be shown, by means of the Schmidt decomposition 



^ I*) = J2a VAS|a{m})|a{«-m}), that 

a<a' a 

(4) 

where represent the eigenvalues of the reduced density 
P{m} = Tr{„_,„}p of the m selected qubits (the same as 
those of P{n-m} when p^ = p). Eq. (g]) differs from the 
entanglement of formation [l^ E[p\ oc — log2 
but is as well a measure of the disorder of the reduced sys- 
tem (see Appendix), satisfying iVj, [p] > Np[p'] if ^ A"', 
where -< denotes "majorized by" (or "more mixed than" ) 

0: K < ELi ^'a for fc = 1, . . . where A^, 

A'^ are sorted in decreasing order and dm = 2™ is the 
subsystem dimension (m < n/2). The formal maximum 
is then obtained for a uniform distribution A^J, — 1/dm, 
in which case Np[p] = {dm — l)/2. Due to convexity, this 
is also the maximum for non-pure states. 

In the nearest neighbor cyclic chain, global negativi- 
ties Np[p] depend on the number m of selected qubits 
and on their spacings. For instance, in a 3-qubit chain 
abc there is single distinct global negativity Na-bc (= 
Nb-ac — Nc~ba), whereas for n = 4 qubits abed, there are 
three distinct negativities Na-bcd, Nab-cd and Nac-bd, 
which measure respectively the entanglement between: 
one qubit and the rest, adjacent pairs, non-adjacent pairs. 
Each negativity Np[p{T)] will have its own limit temper- 
ature Tp above which it vanishes, although the behavior 
for T < Tp may not be monotonous (even entanglement 
vanishing plus reentry may occur for T < Tp jssl l36l|). A 
remarkable feature of these global limit temperatures is 
that for Hamiltonians of the form ffl, they are strictly 
independent of the uniform field b [25| (see Appendix), 
even though A'p[p(T)], as well as the entanglement of the 
ground state, depend of course on b. 

The bipartite entanglement of subsystems of to < n 
qubits can be analyzed in a similar way by evaluation 
of the reduced negativities Npi[p^m}] determined by the 
reduced density p{„i}, with p' = {k} — {m— k} a biparti- 
tion of the subsystem. These negativities will in general 
also depend on the relative location of the k qubits, and 
satisfy inequalities of the form (20| N^ig^_^m-k}[P{m}] ^ 
N{k}~{rn+i-k}[P{7n+i}], as tracing out one qubit is a 
LOGG operation. The associated limit temperatures will 
then satisfy Tik}-{m-k} < T{k}-{m-i-i-k}, being thus 
higher for global bipartitions. In particular, the neg- 
ativity of one qubit with the rest (Na-*) is an upper 
bound to all pairwise negativities Na^b, A^a-c, etc., so 
that Ta-* is an upper bound to all pairwise limit tem- 
peratures. Limit temperatures Tpi of reduced negativities 
^p'[P{m}iT)] depend in general on the field b. 

For instance, for a three-qubit cyclic chain there is a 
single reduced two-qubit density pab = Trcp (any other 
choice of pair is equivalent), with negativity Na-b[pab] < 
Na-bc[p\i while for a four qubit chain there is one dis- 
tinct three-qubit density pabc = Tr^ p, with two different 
negativities Na-bc[Pabc], Nb-ac[Pabc], satisfying Na-bc < 
{Na-bcd, Nab-cd), Nb-ac < {N a- bed, Nac-bd) ■ There are 
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also two different pair densities pab = T^^cdP, Pac = TT^bdP, 
whose negativities Na-t, Na-c, measure the entangle- 
ment of contiguous and non-contiguous qubits and satisfy 

Na-b < {Na-bc,Nb^ac), Na-c < Na-bc- 



III. RESULTS 

We consider in what follows a cyclic chain with near- 
est neighbor coupling, and define the anisotropy as A = 
fz/l^l, with V = Vx. At T = and fixed 6^0, the ground 
state of the Hamiltonian will experience a series of [n/2] 
transitions \M\ — >■ |M| -I- 1 as A increases, starting from 
\M\ = (1/2) for n even (odd) and large negative A, and 
ending in an aligned state with maximum spin \M\ — n/2 
for A > Ac (6). This state is completely separable while 
all ground states with \M\ < n/2 are entangled, so that 
Ac (6) indicates the entangled-separable border at T = 0. 
The energy of the aligned state |^'o) (= \i ... i) for 
b > 0) is Eq ^ -n{2\b\ + v^)/4:, while for v > 0, that 
of the lowest \M\ = 1 state, which is a W^-type state 
(l*i> « Eti s;l*o) for b > 0), is El = Eo + \b\+v,-v, 
so that the transition occurs at 



Ac(fe) = 1 - IV^^I (w > if n odd) , 



(5) 



Eq. ([U is of course also valid for w < if n is even. 

The entanglement borders for T > are, however, 
quite different. As limit temperatures for non-zero global 
negativities are independent of the field b, so will be 
the concomitant borders A^;(r). Global bipartite entan- 
glement for r > and 6 7^ will then also arise for 
Ac(6) < A < 1. Moreover, we will see that it may also 
arise for A > 1 if n > 4. 

Two and three qubit case. Let us first discuss the be- 
havior for two and three qubits, where analytical expres- 
sions for r > can be found (for 71 = 2 we consider just a 
single term v^^ in ([T])). H can be rewritten in these cases 
in terms of the total spin components Sa = X]"=i ■^q 



H = bS,- ^JviSl + S^) + v,Sl] + So , 



(6) 



where Eq = n{2v + Vz)/i, so that its eigenstates have 
good total spin S with energies 

eIj = bM - \[vS{S + M\v, - v)] + Eo , (7) 

where 5 = 0, 1, for n = 2, and S* = 1/2 (two-fold de- 
generate), 3/2, for n — 3, with < S. In the fully 
connected case Eqs. (H))-© are of course valid for all n, 
with S = (1/2), . . . , n/2 for n even (odd). 

For n — 2, the eigenstates of H are then either separa- 
ble (the aligned states {\SM) — |1, ±1)) or maximally en- 
tangled (the BeU states \SM) = |10), |00) cx \U) ± lit)), 
so that the ground state is maximally entangled for 
A < Ac(5) (with = 1 if w > and 5 = if i; < 0). On 
the other hand, any two-qubit mixed state of the form 
p = J2s M Pfi\^'^-^)i^'^^\ entangled whenever 



(see Appendix), which in the thermal case ([2]) (p' 
e~^Ml^ /Z) leads to the 6-independent border 

A < l-2tln[2/(l-e-i/*)], t~T/\v\. 



(9) 



Thermal entanglement arises then for A < 1 V 6, im- 
plying reentry for T > if Ac (6) < A < 1. The limit 
temperature Ta_6(A) determined by (jH)) is a decreasing 
function of A (see Fig. [1]), that vanishes for A 1 and 
diverges for A -00 (A < 1 - 2tln{2t) for t > 1) . 

For n — 3 qubits, the behavior of (7lo6aZ bipartite entan- 
glement is qualitatively similar to that for n = 2 (Fig. 1). 
For V > 0, the ground state is a W^-type entangled state 

if A < Ac(6) m {\SM) = i§,-i) cx im)+i;n> + i;;t) 

for & > 0), with global negativity Na-bc — V2/3 ~ 0.47 
and pair negativity Na-b = (1 — V5)/6 ~ 0.21, becom- 
ing the aligned state ||,±|) for A > Ac (6). Let us add 
that if 6 = 0, the states ||,±i) become degenerate but 
their mixture ^ Em=±i/2 (the T ^ limit 

of p{T) for A < Ac(5)) remains entangled, with lower 
values Na-bc = (VS - l)/3 w 0.24, Na-b = 1/6. 

On the other hand, for any 3-qubit mixed state p = 
J2s mPm-^m^ where P^j denotes the projector onto the 
subspace with total spin S and z-component M, the 
global negativity Na-bc will be non-zero if and only if 
(see Appendix) 



\pI%-pI%\ > \/^pli%(p"^w2+p%/^)^ (10) 

for 1/ = 1 or V = —I, which in the thermal case implies 



A < 1 - tlnl 



3(2- 



-3/2t\ 



2(l-e-3/2*)2J 



(v > 0) 



(11) 



Global entanglement is then again feasible for A < 1 
V & (A < 1 — tin 3 for t ^ 1) and the ensuing limit 
temperature Ta-bc is a decreasing function of A, that 
vanishes for A — > 1 and diverges for A — >■ —00 (A < 
1 - 2t ln(V2t) if t > 1) (see Fig. [J). 

The behavior of the pairwise limit temperature Ta-b 
for 71 = 3 is however quite different [29!, ^33] . The reduced 
two-qubit density pa-b is entangled in a smaller region 
which depends on the field b and is determined by the 
equation (see Appendix) 

A < 1 - tin ■ 



\pl-p'o\>2jp\pli, 



(8) 



^^2(^2 _ 1) + 2(1 + ?7)(1 - a)2 - 777 

(12) 

where 7 = l-|-2a, a — e ^Z^* and 77 — cosh(6/T). For t <C 
1,A<1 — tln3as before, so that for T > 0, pa-b is also 
entangled for A < 1 V 6. However, the denominator in 
(|12p vanishes at a finite temperature tc{b), implying that 
Ta-b/v approaches a finite limit tc{b) for A — >■ —00, in 
contrast with Ta-bc- Hence, for large negative A there is 
a large range of temperatures where only global bipartite 
entanglement persists, with Ta-bc/Ta-b ~^ 00 for A — 
—00. The limit tc{b) is an increasing function of with 
tc(0) = 3/(41n2) sa 1.08. 
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FIG. 1: (Color online). Limit temperatures Ta-bc for global 
bipartite entanglement for n — 3 qubits as a function of the 
asymmetry A — Vz/\v\, for both signs of v (they are inde- 
pendent of the uniform field b). Also depicted are the corre- 
sponding limit temperatures for pairwise entanglement Ta-t 
for 6 = (that for v < vanishes), and the concomitant result 
for n = 2 qubits (the same for both signs of v and any b). 

For V < 0, Eq. PT|) is to be replaced by 

3(1 + 2e-3/2*) 
^ < V2 - nn[X±__i] („<o), (13) 

so that global bipartite entanglement starts for A < 1 /2 
(A < 1/2 - iln(3/2) for t < 1), although for t > 1 
Eq. ([T3|) and (jlip arc almost coincident, implying that 
Ta-bc is almost the same as for w > for large nega- 
tive A (Fig. [T|). This is not the case for Ta-b, which 
for u < is determined by Eq. p2)) with 7 = 2 + a, 
and is lower than the value for u > 0, vanishing for 
& — >■ (no pairwise entanglement in the absence of field 
[29I [13]). This effect is due to the larger degeneracy 
present for v < 0, where the ground state corresponds to 
S = \M\ = 1/2 for A < A- (6) = 1/2 - \v/b\, being then 
two-fold degenerate for 7^ 0, and to an aligned state 
(|M| = 3/2) if A > A~{b). For A < A-{b) and b ^ 0, 
the mixture p = \Pi/2,±i/2 (the T limit of p{T)) 
is still fully entangled, with Na-bc = \/2/6 ~ 0.23 and 
Na-b = (\/2 - l)/6 « 0.07. However, for 6 = 0, states 
\S,ztM) become degenerate and the ensuing T 
limit, p — j{Pi/2.i/2 + Pi/2,-1/2), has still global entan- 
glement (Na-bc — 1/6), but no pairwise entanglement 
{Na-b = 0), explaining the vanishing of Ta-b in this case. 

Four qubit case. A surprise comes already for the n = 4 
chain (Fig. [2]), where for T > 0, global bipartite entangle- 
ment is seen to arise also for A > 1. The ground state is 
entangled just for A < Ac (6) [Eq. ([S])], with the transi- 
tions \M\ \M\ + l occurring at A = {l-2b-¥) / {l+b), 
b = \b/v\ (0 1), and A = Ac(fe) (1 2), collapsing 
both into a single 2 transition at A = 1 for 6 = 0. 
The corresponding limit temperatures are depicted in the 
top panel, where it is seen that they all vanish for A — > 1, 
but those corresponding to Nab-cd (entanglement be- 
tween contiguous Y>'AV[s) and Na-bcd (entanglement of one 
qubit with the rest) become again non-zero for A > 1, in- 
dicating the reentry of the corresponding entanglement. 
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FIG. 2: (Color online). Top: Limit temperatures for global 
{Tab-cd, Tac-bd, Ta-bcd), pairwisc {Ta-b, Ta-c) and three- 
qubit {Ta-bc, Tb-ac) negativities for n = 4 qubits in a near- 
est neighbor cyclic chain, as a function of the asymmetry A. 
There is global entanglement also for A > 1. Results for re- 
duced systems correspond to zero field. Center: The behavior 
of the global negativity Nab-cd as a function of temperature 
and A (with difi'erent scales for A < 1 and A > 1) for 6 = 0. 
Bottom: Negativities of the state H14p as a function of A. For 
6 = 0, they represent the ground state negativities for A < L 

The highest limit temperature for A < 1/2 corresponds 
to Nac-bd (entanglement between non- contiguous pairs), 
but changes to Nab-cd for A > 1/2. The thermal behav- 
ior of Nab-cd for 6 = is depicted in the central panel, 
where it is seen that the reentry for A > 1 is actually 
weak and vanishes smoothly for T — )• 0. Let us remark 
that no entanglement reentry for A > 1 takes place in 
the analogous fully connected four qubit case. 

The M = ground state in the four qubit cyclic chain 
is an entangled state of the form 

\^) = a(itm) + ■•■) + mm) + ■■■), (14) 

where . . . denotes similar states obtained by translation 
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and PI a = (V8 + A2 - A)/2, with ^a^ + 2^^ = 1 and 
/3 > a for A < 1. Its energy is —\v\P/a. This state 
exerts a strong influence on the entanglement of p{T) 
even for A > 1. It has maximum entanglement between 
one qubit with the rest, Na-bcd = 1/2 (the reduced one- 
qubit density pa is maximally mixed), while Nac-bd — 
p{Aa + 13), and N^b-cd = N^c-bd if A < 1 and Ga^ - /J^ 
if A > 1 (Fig. [21 lower panel), becoming Nab-cd much 
larger than Nac-bd for large A. This seems to affect the 
most persistent negativity of p{T) for A > 1, allowing a 
positive value of Nab-cd and reducing Nac-bd to 0. The 
|M| = 1 ground states are VF-states (oc | till) + I itii) + 
. . . for M = -1), with lower values Nab-cd = Nac-bd = 
1/2 and Na-bcd = \/3/4. 

We also depict in the upper panel the limit tempera- 
tures for 6 = of all reduced negativities, namely Ta-bc, 
Tb-ac for the three qubit density pabc and Ta-b, Ta-c for 
the reduced pair densities pab and pac- They all vanish 
for A > 1, so that in this region just global bipartite en- 
tanglement persists (a result valid for any b). It is also 
seen that while Ta-b (adjacent qubits) increases as A 
decreases, diverging for A — — oo (in contrast with the 
behavior for n = 3), Ta-c (non-adjacent qubits) vanishes 
for A < 0. This is a consequence of the ground state 
behavior of Na-c, which vanishes for A < (see bottom 
panel). The state (dH) leads to Na-b = a(2^ - a) for 
A < 7/2 (and for A > 7/2), and Na-c = ^a^ - 
for A > (and for A < 0), with Na-b > Na-c for 
A < 1. The ordering of negativities is of course in agree- 
ment with the discussion of sec. II. 

The behavior of global negativities for n = 5 qubits 
and V > (Fig. [31 upper panel) is quite similar, although 
there are some important changes: a) For T > 0, all 
global negativities exhibit a reentry for A > 1, including 
that of the non-contiguous 2 — 3 partition, Nac-bde, with 
Tac-bde lying close to Tab-cde (contiguous 2-3 partition) 
for A > 1/2; b) The negativity of one-qubit with the 
rest, Na-bcde, remains non-zero V A if T > 0, includ- 
ing the isotropic case A = 1, providing the highest limit 
temperature for 0.85 A < 1.35; c) For A > 1 there 
is no pairwise nor three qubit entanglement at any field, 
but there is four qubit entanglement (all reduced four- 
qubit negativities are small but non-zero). Besides, for 
A — > — oo the limit temperature for adjacent pairwise en- 
tanglement Ta-b approaches a finite limit [23|, as occurs 
for n = 3, exhibiting a maximum at A ss —5.6 for 6 = 0, 
but all those for three and four adjacent qubits {Ta-bc, 
Tb-ac, Ta-bcd, ■ ■ •) diverge, as the global ones. Note also 
that Ta-c vanishes below a certain limit (A < 0.1 for 
6 = 0), reflecting the vanishing of Na-c at T = below 
this value. For n = 5 and > 0, the ground state transi- 
tions |M| \M\ + 1 are located at A = (l-6-62)/(l-|-6) 
(i — |) and Ac (6) (| — )• |), collapsing into a single 
i — >■ I transition at A = 1 if 6 = 0. 

Global limit temperatures for rt = 5 and v < (lower 
panel) are quite close to the previous ones. The main 
difference is that all of them remain non-zero for all A, 
the most persistent negativity corresponding to N ac-bde 
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FIG. 3: (Color online). The set of limit temperatures for 
global negativities, Ta-bcde, Tab-cde, Tac-bde, for 71 = 5 qubits 
in the cyclic chain as a function of A for v > (top) and v < 
(bottom). The concomitant limit temperatures for pairwise 
{Ta-b, Ta-c) and the most persistent three-qubit {Tb-ac) neg- 
ativity for 6 = are also depicted. Ta-c vanishes for i; < 
and 6 = 0. 

if A < 0.5 and Nab-cde if A > 0.5. Limit temper- 
atures for pairwise entanglement exhibit more signifi- 
cant differences. For 6 = there is no non-contiguous 
pairwise entanglement if w < 29] {Ta-c ~ 0), since 
the ground state is four-fold degenerate and leads to 
Na-c = V A, while Ta-b is lower (saturating again 
for A — > — oo). For v < 0, the ground state transitions 

\M\ ^ \M\ + 1 occur at A = -i^^iZgM±^^ and 

A = A^(6) = j{l + VS) - 6, with both coUapsing at 
A-(0) « 0.81 for 6 = 0. Again, for A > Ac(0) there is 
no two- nor threc-qubit negativity, but there is four qubit 
entanglement. 

In order to appreciate the trend for larger n, Fig. 
[H depicts the set of limit temperatures of the 17 dis- 
tinct global negativities existing for n = 8 qubits 
ab . . . gh, consisting of Na-* (one qubit with the rest), 
and those between different pairs and the rest {Nab-*, 
Nac-*, Nad-*, Nae-*), three qubits and the rest {Nabc-*, 
Nabd-*, Nabc-*, Nace-*, Nacf-*) and four qubits and 

the rest {Nabcd-*, Nabce-*, Nabcf-* = Nabde-*, Nabef-*, 

Nabdf-*, Nabdg-*, Naceg-*)- All thcsc tcmpcratures re- 
main non-zero VA, including the region A > 1 where the 
ground state is aligned for any 6, confirming the tendency 
observed for n = 4 and 5. The highest limit tempera- 
ture corresponds to the full non-contiguous 4-4 partition 
{Taceg-*) for A < 0.5, but changes to that between non- 
contiguous pairs of adjacent qubits {Tabcf-*) for A > 0.5 
(in this region Taceg-* becomes the lowest global limit 
temperature). Again, the "depth" of the entanglement 
for A > 1 is up to four qubit subsystems, being zero all 
pairwise and three-qubit negativities. The limit temper- 
ature for adjacent pairwise entanglement for 6 = lies 
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FIG. 4: (Color online). The limit temperatures of all seven- 
teen global negativities for n = 8 qubits in a cyclic chain as a 
function of asymmetry. The highest and lowest temperatures 
in each sector are specially indicated. The limit temperatures 
for adjacent and non-adjacent pairwise entanglement for 6 = 
are also depicted. 

well below the bundle of global limit temperatures, while 
those of non-adjacent pairs are small and non-zero just 
in a small interval between A = and A = 1, where the 
ground state negativity is non-zero. Let us also remark 
that the revival for A > 1 is sensitive to the interaction 
range. For instance, there is no entanglement revival in 
the analogous fully connected case. 

Finally, Fig. ([5]) compares results for the entanglement 
between one qubit with the rest for different n, assum- 
ing V > 0. As previously stated, Na-* and hence Ta-* 
are upper bounds to all pairwise negativities and limit 
temperatures. It is first seen in the upper panel that 
the global limit temperatures T^-* for different n rapidly 
converge to a common value, results for n = 6 — 10 be- 
ing undistinguishable in the scale and range of the figure. 
This temperature exhibits a slope discontinuity around 
A « 0.8, where the number k of negative eigenvalues of 
the partial transpose exhibits a minimum (while in the 
case of pure states fc = 1 for Na-* [Eq. ([4])], it can be 
much larger than 1 for non-pure states). 

The corresponding global negativities A^a-* also ap- 
proach a limit curve for A < 1, as seen in the central 
panel for 6 = and A = — 1. Note that for all even n, 
Na-* starts at its maximum value 1/2 at T = 0, since 
the reduced one-qubit density pa for the M = ground 
state is maximally mixed and Eq. ^ reaches then its 
maximum value. At fixed T > 0, they tend initially to 
decrease for increasing n, rapidly approaching the limit 
curve. On the other hand, for n odd and 6 = 0, the 
ground state for A = — 1 has M — ±1/2 and is degener- 
ate, so that the initial value is lower than in the even case 
(see the comment for n — 3), although it increases as n 
increases and approaches the limit curve (for 6^0, Na-* 
would start at ~ 0.47 for n = 3, increasing then with n). 
The lower panel depicts the reentry of the negativity for 
T > at A = 1.5 for increasing n. The behavior is less 




-1 A 1 2 




0.1 0.2 0.3 



T/lvl 

FIG. 5: (Color online) Top: Limit temperatures for non-zero 
negativity between one qubit and the remaining n — 1 qubits, 
for V > and n — 2, ... 10. Center and bottom: The thermal 
behavior of these negativities for A = — 1 and 1.5, for fe = 0. 

monotonous than in the previous case, although the dif- 
ference between results for neighboring odd-even values 
of n tend to decrease as n increases. 



IV. CONCLUSIONS 

The present results demonstrate that limit tempera- 
tures for global bipartite entanglement in XXZ chains 
may differ considerably from those limiting pairwise en- 
tanglement between two qubits in the same chain. In 
contrast with the latter, global limit temperatures are in- 
dependent of the field 6 and do not saturate for A — > — oo 
in odd chains. Moreover, global thermal entanglement in 
nearest-neighbor cyclic chains may naturally arise at low 
levels even for A > 1 if n > 4 (and A > 1 if n > 5), 
i.e., in anisotropic and isotropic chains with fully sep- 
arable ground states, generating entangled bipartitions 
of the whole system as well as of subsystems with no 
less than four qubits. The ordering of global limit tem- 
peratures also changes as A increases from negative to 
large positive values. Many of our results rapidly sat- 
urate as n increases (Fig. [5]), indicating that they may 
remain stable for chains with a larger number of qubits. 
The present study provides thus a more comprehensive 
understanding of the borders of thermal entanglement 
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in XX Z spin chains, suggesting that quantum schemes 
based on entanglement between multi-spin parties may 
be more resistant to the effects of thermal randomness. 

Acknowledgment. NC and RR acknowledge support, 
respectively, from CONICET and CIC of Argentina. 

Appendix 

Negativity of pure states. Eq. (|4]) can be rewritten as 

^p[p] = i{[5/(P{™}) + i]'-i}, 

where 

'5'/(P{m}) = Tr/(p{„}) , f[p] = ^-p, 

is a simple non-additive entropy of p^m} fisj - since the 
function / is concave in [0, 1] and satisfies /(O) = /(I) = 
0. As a consequence, > 0, being maximum for 

P{m} = I /dm and zero when p^m} is pure, and satisfies 
Sf{,P{m}) > Sf{p[^}) if P{m} ^ P{,„}, where the last con- 
dition means that the eigenvalues of P{m} are majorized 
by those of p'^^} E3- These properties are then also 
satisfied by the pure state negativity Np since it is an 
increasing function of S*/ that vanishes for 5/ = 0, and 
ensure that it cannot increase under LOCC. 

Independence of global limit temperatures from the uni- 
form magnetic field b. Let us consider a Hamiltonian 

H = bS,+V, 

where V is independent of b and satisfies [V^, S'z] = 0. For 
an arbitrary bipartition p = {m} — {n — m} of an n-qubit 
system and a standard basis of states |M{„j}, Af — M{„}), 
where M is the eigenvalue of the total spin component Sz 
and M{m} that for the first subsystem (remaining labels 
omitted), the partial transpose of the exponential D{b) = 
exp[-iJ/r] = exp[-bSz/T]exp[-V/T] in the previous 
basis satisfies [l^ 

D*''{b) = exp[^bSz/2T]D'''{0) exp[~bSz/2T], 

where _D*p(0) has matrix elements between states with 
total spin component M and M ± 2k, with k integer. 
While the negativity will depend on b, the limits for the 
positivity of D*p{b) will not, since exp[—bSz/2T] is real 
and positive, being the same as those for Z?*p(0). Global 
limit temperatures will then be independent of b. This 
applies in particular to any XXZ type Hamiltonian. 

Negativity of mixtures of two- and three-qubit states 
with good angular momentum. The negativity iVa-6 of 
any two qubit mixed state of the form 

s 

P= E E Pii\SM){SM\, (15) 

S=0,1 M=-S 



where | SM) denote the states with good total spin S and 
component M, can be shown to be [3a] 

Na-b = iMax[^(pi-pO)2 + (pi-pLJ2 _pl_ pl^^Q^ 

(16) 

The condition Na-t > leads then to Eq. ([8]). Note 
also that the single qubit density is pa — Trbp = 
J2u=±l/2l'^W){^\^ with =pI + i(pj] +pj). 

For n — 5 qubits, a set of 8 orthonormal states with 
good total spin 5* and spin component M is given by 
III) = Ittt), IM) = (litt) + int) + IUt))/V3, |iia) = 

(iin)-int))/V2, jiifo) = (i;n)+int)-2im))/V6, 

and the corresponding partners for M — —M. For a 
three-qubit density of the form 

P= E E PmPm, (17) 

5=1/2,3/2 M=-S 

where denotes the projector onto the subspace with 
total spin S and spin component M, an analytic eval- 
uation of the eigenvalues of the partial transpose of p 
yields the following expression for the global negativity 

Na-bc — Nab-c = Nca-b, 

Na-bc = ^ E ^^^^[ 
//o 3/2 o 3/2 1/2n2 i o/ 3/2 1/2 x, 

-ip%l2-'^P%-P%)A- (18) 

The condition Na-bc > leads then to Eq. The 
ensuing reduced two qubit density Pa-b = Trcp (identical 
to Pa-c, Pb-c) has again the form with 

1 /n 3/2 ,3/2 , r, 1/2 N /o 

P±l = (3p±'3/2+P±i/2 + 2p±i/2)/3, 

pI = i2{p'/;,+p%)+p',)/3, p"o=pI^/1+P% 

The condition ([5]) for pairwise entanglement leads then 
to 

I E p%i-pj;2\ > n m%+p%i+'pI%- (19) 

u=±l u=±l 

In the thermal case for the Hamiltonian ([1]), Eq. 
leads to Eq. 
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